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ABSTRACT 

Using a large galaxy group catalogue constructed from the Sloan Digital Sky Survey Data Release 
4 (SDSS DR4) with an adaptive halo-based group finder, we investigate the luminosity and stellar 
mass functions for different populations of galaxies (central versus satellite; red versus blue; and 
galaxies in groups of different masses) and for groups themselves. The conditional stellar mass function 
(CSMF), which describes the stellar distribution of galaxies in halos of a given mass for central and 
satellite galaxies can be well modeled with a log-normal distribution and a modified Schechter form, 
respectively. On average, there are about 3 times as many central galaxies as satellites. Among the 
satellite population, there are in general more red galaxies than blue ones. For the central population, 
the luminosity function is dominated by red galaxies at the massive end, and by blue galaxies at the 
low mass end. At the very low-mass end (M» < 10 9 h~ 2 M©), however, there is a marked increase in 
the number of red centrals. We speculate that these galaxies are located close to large halos so that 
their star formation is truncated by the large-scale environments. The stellar-mass function of galaxy 
groups is well described by a double power law, with a characteristic stellar mass at ~ 4 x 10 10 ft~ 2 M©. 
Finally, we use the observed stellar mass function of central galaxies to constrain the stellar mass - 
halo mass relation for low mass halos, and obtain M* iC oc M^' 9 for Mh <§C 10 11 /i _1 M©. 
Subject headings: dark matter - large-scale structure of the universe - galaxies: halos - methods: 
statistical 



1. INTRODUCTION 

In recent years, great progress has been made in our 
understanding about how galaxies form and evolve in 
dark matter halos owing to the development of halo mod- 
els and the related halo occupation models. The halo oc- 
cupation distribution (hereafter HOD), P(N\Mh), which 
gives the probability of finding N galaxies (with some 
specified properties) in a halo of mass Mh, has been ex- 
tensively used to study the galaxy distribution in dark 
matter halos and galaxy clustering on large scales (e.g. 
Jing, Mo & Borner 1998; Peacock & Smith 2000; Sel- 
jak 2000; Scoccimarro et al. 2001; Jing, Borner & Suto 
2002; Berlind & Weinberg 2002; Bullock, Wechsler & 
Somerville 2002; Scranton 2002; Zehavi et al. 2004, 2005; 
Zheng et al. 2005; Tinker et al. 2005; Skibba et al. 2007; 
Brown et al. 2008) . The conditional luminosity function 
(hereafter CLF), $(£|M h )d£, which refines the HOD 
statistic by considering the average number of galax- 
ies with luminosity L ± dL/2 that reside in a halo of 
mass Mh, has also been investigated extensively (Yang, 
Mo & van den Bosch 2003; van den Bosch, Yang & Mo 
2003; Vale & Ostriker 2004, 2008; Cooray 2006; van den 
Bosch et al. 2007) and has been applied to various red- 
shift surveys, such as the 2-degree Field Galaxy Redshift 
Survey (2dFGRS), the Sloan Digital Sky Survey (SDSS) 
and DEEP2 (e.g. Yan, Madgwick & White 2003; Yang 
et al. 2004; Mo et al. 2004; Wang et al. 2004; Zehavi 

1 Shanghai Astronomical Observatory, the Partner Group of 
MPA, Nandan Road 80, Shanghai 200030, China; E-mail: 



xhyang@shao.ac.cn 
2 Department of 



University of Massachusetts, 



Astronomy, 
Amherst MA 01003-9305 

3 Max-Planck-Institute for Astronomy, Konigstuhl 17, D-69117 
Heidelberg, Germany 

4 Joint Institute for Galaxy and Cosmology (JOINGC) of Shang- 
hai Astronomical Observatory and University of Science and Tech- 
nology of China 



et al. 2005; Yan, White & Coil 2004). These investi- 
gations demonstrate that the HOD/CLF statistics are 
very powerful tools to establish and describe the con- 
nection between galaxies and dark matter halos, provid- 
ing important constraints on various physical processes 
that govern the formation and evolution of galaxies, such 
as gravitational instability, gas cooling, star formation, 
merging, tidal stripping and heating, and a variety of 
feedback processes, and how their efficiencies scale with 
halo mass. Furthermore, they also indicate that the 
galaxy/dark halo connection can provide important con- 
straints on cosmology (e.g., van den Bosch, Mo & Yang 
2003; Zheng & Weinberg 2007). 

However, as pointed out in Yang et al. (2005c; here- 
after Y05c), a shortcoming of the HOD/CLF models is 
that the results are not completely model independent. 
Typically, assumptions have to be made regarding the 
functional form of either P(N\M h ) or <&{L\M h ). More- 
over, in all HOD/CLF studies to date, the occupation 
distributions have been determined in an indirect way: 
the free parameters of the assumed functional form are 
constrained using statistical data on the abundance and 
clustering properties of the galaxy population. An al- 
ternative method that can directly probe the galaxy - 
dark halo connection (e.g. HOD/CLF models) is to use 
galaxy groups as a representation of dark matter halos 
and to study how the galaxy population changes with 
the properties of the groups (e.g., Y05c; Zandivarez et 
al. 2006; Robotham et al. 2006; Hansen et al. 2007; 
Yang et al. 2008). For such a purpose, one has to prop- 
erly find the galaxy groups that are closely connected 
to the dark matter halos. In recent studies, Yang et al. 
(2005a; 2007) developed an adaptive halo-based group 
finder that has such features 5 . This group finder has 

5 In this paper, we refer to systems of galaxies as groups regard- 
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been applied to the 2dFGRS (Yang et al. 2005a) and to 
the SDSS (Weinmann et al. 2006a; Yang et al. 2007). 
Detailed tests with mock galaxy catalogues have shown 
that this group finder is very successful in associating 
galaxies according to their common dark matter halos. 
In particular, the group finder performs reliably not only 
for rich systems, but also for poor systems, including iso- 
lated central galaxies in low mass halos. This makes it 
possible to study the galaxy-halo connection for systems 
covering a large dynamic range in masses. With a well- 
defined galaxy group catalogue, one can then not only 
study the properties of galaxies in different groups (e.g. 
Y05c; Yang et al. 2005d; Collister & Lahav 2005; van 
den Bosch et al. 2005; Robotham et al. 2006; Zandivarez 
et al. 2006; Weinmann et al. 2006a, b; van den Bosch et 
al. 2008; Mcintosh et al. 2007; Yang et al. 2008), but 
also probe how dark matter halos trace the large-scale 
structure of the Universe (e.g. Yang et al. 2005b, 2006; 
Coil et al. 2006; Berlind et al. 2007; Wang et al. 2008a). 

Recently, this group finder has been applied to the 
Sloan Digital Sky Survey Data Release 4 (SDSS DR4), 
and the group catalogues constructed are described in de- 
tail in Yang et al. (2007; Paper I hereafter). In these cat- 
alogues various observational selection effects are taken 
into account, and each of the groups is assigned a re- 
liable halo mass. The group catalogues including the 
membership of the groups are available at these links 6 
'. In Yang et al. (2008; Paper II hereafter) we have 
used these group catalogues to obtain various halo occu- 
pation statistics and to measure the CLFs for different 
populations of galaxies. In this paper, the third in the 
series, we will focus on the conditional stellar mass func- 
tions (CSMFs) for different populations of galaxies. In 
addition, we will also examine the general luminosity and 
stellar mass functions for different populations of galaxies 
and for groups themselves. Finally, we will demonstrate 
how to use the observed luminosity and stellar mass func- 
tions for central galaxies to constrain the HOD in small 
halos. 

This paper is organized as follows: In Section [2] we 
describe the data (galaxy and group catalogues) used in 
this paper. Section [4] presents our measurement of the 
CSMFs for all, red and blue galaxies. Sections [3] and [5] 
present our measurement of the luminosity and stellar 
mass functions for galaxies and groups, respectively. In 
Section [6] we probe the properties of the central galax- 
ies that can be formed in those small halos. Finally, 
we summarize our results in Section [7] Throughout this 
paper, we use a ACDM 'concordance' cosmology whose 
parameters are consistent with the three-year data re- 
lease of the WMAP mission: Q, m = 0.238, fl A = 0.762, 
n s = 0.951, h = 0.73 and a 8 = 0.75 (Spergel et al. 2007). 
If not quoted, the units of luminosity, stellar and halo 
masses are in terms of h~ 2 L Q , h~ 2 M Q and M Q , re- 
spectively. Finally, unless noted differently, the luminos- 
ity functions and stellar mass functions arc presented in 
units of /i 3 Mpc~ 3 dlogL and /i 3 Mpc~ 3 d log M* , respec- 
tively, where log is the 10 based logarithm. 

less of their richness, including isolated galaxies (i.e., systems with 
a single member) and rich clusters of galaxies. 

6 http:/ /gax. shao.ac.cn/data/Group. html 

7 http:/ /www. astro. umass.edu/~xhyang/Group. html 



2. DATA 
2.1. Galaxy and group catalogues 

The data used in our analysis here are the same as 
those used in Paper II. Readers who have already read 
through Paper II may go directly to the next subsection. 

The group catalogues are constructed from the New 
York University Value-Added Galaxy Catalogue (NYU- 
VAGC; Blanton et al. 2005b), which is based on the 
SDSS Data Release 4 (Adelman-McCarthy et al. 2006), 
but with an independent set of significantly improved 
reductions. From NYU-VAGC we select all galaxies in 
the Main Galaxy Sample with redshifts in the range 
0.01 < z < 0.20 and with a redshift completeness 
C > 0.7. As described in Paper I, three group samples 
are constructed from the corresponding galaxy samples: 
Sample I, which only uses the 362356 galaxies with mea- 
sured r-band magnitudes and redshifts from the SDSS, 
Sample II which also includes 7091 galaxies with SDSS 
r-band magnitudes but redshifts taken from alternative 
surveys, and Sample III which includes an additional 
38672 galaxies that lack redshifts due to fiber collisions 
but that are assigned the redshifts of their nearest neigh- 
bors. Although this fiber collision correction works well 
in roughly 60 percent of all cases, the remaining 40 per- 
cent are assigned redshifts that can be very different from 
their true values (Zehavi et al. 2002). Samples II and III 
should therefore be considered as two extremes as far 
as a treatment of fiber-collisions is concerned. Unless 
stated otherwise, our results are based on Sample II. For 
comparison, we also present some results obtained from 
Sample III. 

The magnitudes and colors of all galaxies are based on 
the standard SDSS Petrosian technique (Petrosian 1976; 
Strauss et al. 2002), have been corrected for galactic ex- 
tinction (Schlegcl, Finkbeiner & Davis 1998), and have 
been i^-corrected and evolution corrected to z = 0.1, us- 
ing the method described in Blanton et al. (2003a; b). 
We use the notation 01 M r — 51og/i to indicate the re- 
sulting absolute magnitude in the r-band. The galaxies 
are separated into red and blue subsamples according to 
their bi- normal distribution in the (g — r) color (Baldry 
et al. 2004; Blanton et al. 2005a; Li et al. 2006), using 
the separation criteria (see Paper II), 



o.i 



= 1.022 - 0.0651x - 0.00311^ 



(1) 



where x= 01 M r - 5 \ogh + 23.0. 

Stellar masses, indicated by M*, for all galaxies are 
computed using the relations between stellar mass-to- 
light ratio and (5 — r) color from Bell et al. (2003), 

log fe _ 2 ^ = -0.306+ 1.097 [°-°(.9-r)] -0.10 

-0.4( M,. - 5 log /i- 4.64). (2) 

Here °'°(g — r) and M r — 51og/i are the (g — r) color 
and r-band magnitude K + E corrected to z = 0.0; 4.64 
is the r-band magnitude of the Sun in the AB system 
(Blanton & Roweis 2007); and the —0.10 term effectively 
implies that we adopt a Kroupa (2001) IMF (Borch et 
al. 2006). 

For each group in our catalogue we have two estimates 
of its dark matter halo mass M^: (1) Ml, which is 
based on the ranking of the characteristic group lumi- 
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nosity £19.5 , and (2) Ms, which is based on the rank- 
ing of the characteristic group stellar mass M s t iiar, re- 
spectively 8 . The halo mass is estimated for each group 
with at least one member galaxy that is brighter than 
°- 1 M r — 51og/i = —19.5. As shown in Paper I, these two 
halo masses agree reasonably well with each other, with 
scatter that decreases from ~ 0.1 dex at the low-mass 
end to ~ 0.05 dex at the massive end. Detailed tests us- 
ing mock galaxy redshift surveys have demonstrated that 
the group masses thus estimated can recover the true 
halo masses with a 1-a deviation of ~ 0.3 dex, and are 
more reliable than those based on the velocity dispersion 
of group members (Y05c; Weinmann et al. 2006; Berlind 
et al. 2006; Paper I). Note also that survey edge effects 
have been taken into account in our group catalogue: 
groups that suffer severely from edge effects (about 1.6% 
of the total) have been removed from the catalogue. In 
most cases, we take the most massive galaxy (in terms 
of stellar mass) in a group as the central galaxy (MCG) 
and all others as satellite galaxies. In addition, we also 
considered a case in which the brightest galaxy in the 
group is considered as the central galaxy (BCG). Tests 
have shown that for most of what follows, these two defi- 
nitions yield indistinguishable results. Whenever the two 
definitions lead to significant differences, we present re- 
sults for both. Throughout this paper, results are calcu- 
lated for both samples II and III using both halo masses, 
Ml and Ms- Any significant differences in the results 
due to the use of different samples and mass estimates 
are discussed. 

Finally, we caution that the SDSS pipeline may have 
underestimated the luminosities for bright galaxies (e.g. 
von der Linden et al. 2007; Guo et al. 2009). Accord- 
ing to Guo et al. the NYU-VAGC magnitude is over- 
estimated by about 0.5 ±0.1 at apparent magnitudes 
r - 13.0 and about 0.1 ± 0.1 at r ~ 17.0. Although 
this will not change the halo masses estimated using the 
abundance match to halo mass function, the luminosity 
and stellar mass functions for galaxies and groups are 
shifted slightly at the bright ends if a correction is made 
to the SDSS luminosities. 

2.2. Galaxy and group completeness limits 

Because of the survey magnitude limit, only bright 
galaxies can be observed. This consequently will induce 
incompleteness in the distribution of galaxies with re- 
spect to absolute magnitude and stellar mass, and in the 
distribution of groups with halo mass. In this subsection, 
we discuss such completeness and how to make correc- 
tions. 

As discussed in detail in the Appendix of van den Bosch 
et al. (2008), the apparent magnitude limit of the galaxy 
sample (m r = 17.77) can be translated into a redshift- 
dependent absolute magnitude limit given by 
01 M rilim -51og/i = 

17.77 - DM(z) - fco.i(z) + 1.62(2; - 0.1) . (3) 
where ko.i{z) is the i\~-correction to z = 0.1, the 
1.62(2 — 0.1) term is the evolution correction of Blan- 
ton et al. (2003b), and 

DM(z) = 51ogL> L (z) + 25 (4) 

8 L19.5 and M sto n ar are, respectively, the total luminosity and 
total stellar mass of all group members with 01 M r — 51og/i < 
-19.5. 
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Fig. 1. — The group completeness limit as a function of halo 
mass. Different lines (except the solid one) correspond to different 
combinations of Samples and halo masses estimated (M^ or Mg) 
as indicated. The solid line illustrates a conservative halo mass 
limit for all those combinations. That is the halos below given 
redshift with halo masses larger than M^ij m are complete. 

is the distance module corresponding to redshift z, with 
Dl(z) the luminosity distance in /i _1 Mpc. Using the 
^-corrections of Blanton et al. (2003a; see also Blanton 
& Rowcis 2007), the redshift-dependence is reasonably 
well described by 

k . 1 (z) = 2.51og(^±^y (5) 

At z = 0.1, where (by definition) k .i = — 2.51og(l.l) ~ 
—0.1 for all galaxies (e.g., Blanton & Roweis 2007), this 
exactly gives the absolute magnitude limit of the sample. 
At lower and higher rcdshifts, however, a small fraction 
of the sample galaxies fall below this limit. This owes 
to the fact that fco.i(z) not only depends on redshift but 
also on color. In order to ensure completeness, we use a 
conservative absolute-magnitude limit: 

ai M; jim - 51og/i = ai M r , lim - 5 log h- 0.1 , (6) 

where the term 0.1 takes into account the scatter in the K 
correction. For the stellar masses of the SDSS galaxies, 
we adopt, for given redshift z, the following completeness 
limit: 

log[M* tlim /(h- 2 M Q )] = (7) 

4.852 + 2.246 \ogD L {z) + 1. 123 log(l + z) - 1.186z 

1 - 0.067z ' 

(see van den Bosch et al. 2008) 

Next we consider incompleteness in the actual group 
catalogue. As illustrated in Fig. 6 of Yang et al. (2007), 
the groups within a certain luminosity L19.5 or stellar 
mass Mgtciiar bin (which corresponds to a certain halo 
mass bin) are complete only to a certain redshift. Be- 
yond this redshift the number density of the groups drops 
dramatically. We therefore need to take this complete- 
ness into account. Here we proceed as follows to obtain 
the halo-mass completeness limit at any given redshift. 
First, for a given halo mass Mh, we measure the number 
densities of groups with halo mass within log Mh ± 0.05, 
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n(0.01, Zmax) and n(z max — Az, z max ), within the redshift 
ranges, [0.01, z max ] and [z max - Az, z max ], respectively. 
Here we set Az = 0.005. If the group sample is not 
complete at redshift z max , then the group number den- 
sity n(z max — Az, 2 max ) is expected to drop significantly. 
Starting from z max = 0.2, we iteratively decrease z max 
according to z max = z max — Az, and find the largest red- 
shift z max that satisfies, 

n(z max 

Az, 

-'-max 

) +3ct„ > n(0.01,z max ), (8) 

where a n is the variance in the number density of groups 
within [z max — Az, z max ] among 200 bootstrap samples. 
The z max thus obtained is the one at which the halos are 
complete down to the corresponding halo mass M^. For 
the group samples used in this paper, Samples II and 
III, and for the two sets of halo masses estimated for our 
groups, Ml and Ms, we obtain the value of z max as a 
function of halo mass. The results are shown in Fig. [TJ 
From this plot, we obtain a conservative halo-mass limit, 

logM^un, = {z- 0.085)/0.069+ 12 , (9) 

which is shown as the solid line in Fig. [JJ Clearly, this 
criterion works for both Samples II and III, as well as for 
both Ml and Ms- Thus, for a given redshift z, groups 
with masses > Mh,nm are complete. 

3. THE GALAXY LUMINOSITY AND STELLAR 
MASS FUNCTIONS: CENTRAL VS SATELLITE 

In this section we estimate the luminosity and stellar 
mass functions for different populations of galaxies. Both 
functions have been extensively investigated in the liter- 
ature for galaxies of different colors and morphological 
types (e.g., Lin et al. 1996; Norberg et al. 2002; Madg- 
wick et al. 2002; Blanton et al. 2003b; Bell et al. 2003; 
Fontana et al. 2006). Note that in a very recent pa- 
per, based on the group catalogue constructed from the 
2 degree Field Galaxy Redshift Survey (2dFGRS; Col- 
less et al, 2001) by Tago et al. (2006), Tempel et al. 
(2008) measured luminosity functions for various con- 
tents of galaxies (e.g., central, second ranked, satellite, 
isolated) in groups, as well as for groups. Here we focus 
on the difference between central and satellite galaxies. 
We use direct counting to estimate the luminosity func- 
tion. For each galaxy at a given redshift we calculate its 
absolute magnitude limit according to Eq. Q . If the ab- 
solute magnitude of the galaxy is fainter than this limit, 
it is removed from the counting list. If the galaxy is 
not removed, we first calculate the maximum redshift 
at which the galaxy (with its absolute magnitude) can 
be observed. We then calculate the comoving volume, 
Vcom, between this maximum redshift and a minimum 
redshift z = 0.01. In the counting, the galaxy is assigned 
a weight, 

Wi = 77~~~n ( 10 ) 

where C is the redshift completeness factor in the NYU- 
VAGC. We calculate the luminosity functions for all, red, 
and blue galaxies, respectively. The corresponding re- 
sults are shown in the upper-left panel of Fig. [2j where 
open circles, squares and triangles are the results for all, 
red and blue galaxies, respectively. For clarity the results 
for red and blue galaxies are shifted downwards by a fac- 
tor of 10. By comparing the LFs for red and blue galax- 
ies, one sees that there are more (fewer) red galaxies than 



blue ones at luminosities logL > 9.8 (logL < 9.8). For 
very faint galaxies with logL < 9.0, the red population 
increase dramatically, exceeding that of blue galaxies at 
logL ~ 8.0. 

We can further separate the galaxies into centrals and 
satellites according to their memberships in the groups. 
The corresponding luminosity functions for all, red and 
blue galaxies are shown in the middle-left and lower-left 
panels of Fig. [5J respectively. The color dependence of 
the luminosity function for the centrals resembles that of 
the overall population. However, for the satellite popu- 
lation, the color dependence is somewhat different, espe- 
cially around L = 10 9 h~ 2 LQ, where the luminosity func- 
tion is not suppressed relative to that of blue satellites, as 
is the case for the centrals. One interesting feature in the 
luminosity function of red central galaxies is that there 
are many very faint red central galaxies (slightly more 
than the blue ones) with luminosity logL ~ 8.0. Such a 
population is not expected in the standard galaxy forma- 
tion models, where very small halos are expected to host 
only blue centrals. However, Wang, Mo & Jing (2007) 
found that the large-scale tidal field may effectively trun- 
cate the mass accretion into small halos. If gas accretion 
is also truncated in this process, the central galaxies in 
these small halos are expected to be red. More recently, 
Ludlow et al. (2008) found that some sub-halos that 
have at some time been within the virial radius of their 
main progenitors can be ejected, so that some low-mass 
halos at z = outside/near the larger virialized halos 
may have experienced tidal and ram-pressure stripping, 
so that they have stopped forming stars. In both scenar- 
ios, a population of faint red galaxies is expected to be 
present in the vicinity of high density regions but outside 
large, virialized halos. This population may be responsi- 
ble for the upturn of the luminosity function of red cen- 
tral galaxies at the faint end. In a separate paper (Wang 
et al. 2008b), we will discuss the properties and spatial 
clustering of the faint red population in more detail, and 
check the contamination due to false groups near massive 
ones using mock galaxy and group catalogs. 

We also measure the stellar mass functions separately 
for all, red, blue, central and satellite galaxies, taking into 
account the stellar mass completeness limit described by 
Eq. . Only galaxies with stellar masses above the com- 
pleteness limit are used in the estimate. In the right pan- 
els of Fig. [2j we show the stellar mass functions using the 
same symbols as in the left panels. The general behaviors 
of the stellar mass functions are very similar to the cor- 
responding luminosity functions. For reference, we list 
the luminosity functions in Table [TJ and the stellar mass 
functions in Tabled 

We use the following Schechter function to fit the lu- 
minosity function: 



$(L) = (f>* 
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For the stellar mass function luminosity functions, we use 
a similar model: 



(12) 



For each model, there are three free parameters, the am- 
plitude (f>* , the faint end slope a and the characteristic 
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The galaxy luminosity functions 
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.2813 


3 


1312 


± 


o 


.2737 


3 


4498 


± 


o 


5347 


1 


.0388 


± 


o 


,1995 


2 


,4110 


± 





.4263 


1 


5326 


± 





5281 





8124 


± 


o 


.2851 


0, 


7202 


± 


o 


3201 


8.4 


1 


.5373 


± 





.3176 


1 


.4754 


± 


o 


.1732 


3 


.0620 


± 


o 


2409 


3 


.0877 


± 


o 


4414 





.7873 


± 


o 


1432 


2 


,3004 


± 





,3883 


1 


4496 


± 


o 


,4447 





.6881 


± 


o 


1755 


o 


,7615 


± 





3339 


8.5 


4 


.7646 


± 





.2712 


1 


.4727 


± 


o 


.1820 


3 


2919 


± 


o 


.1862 


3 


.2717 


± 


o 


3598 





8272 


± 


o 


.0833 


2 


,4445 


± 





,3416 


1 


4930 


± 





,3578 





.6455 


± 


o 


1481 


0, 


,8474 


± 


o 


,2594 


8.6 


5 


0611 


± 


0. 


1803 


1 


3903 


± 


o 


.1616 


3 


6708 


± 





2012 


3 


4916 


± 


o 


.3415 


o 


.6852 


± 


o 


0744 


2 


.8064 


± 


o 


3192 


1 


5695 


± 





3034 





7051 


± 


o 


.1422 


0, 


8644 


± 


o 


.2057 


8.7 


4. 


9694 


± 


0. 


.1665 


1 


.5349 


± 





.1161 


3 


.4344 


± 





1741 


3 


.1904 


± 





.3180 





.6093 


± 





.0819 


2 


,5811 


± 


0, 


,2763 


1 


7790 


± 





,2697 





,9256 


± 





.1269 





8534 


± 





1816 


8.8 


4. 


.7649 


± 


0. 


.1504 


1 


.3059 


± 





.1176 


3 


.4589 


± 


0. 


1882 


3 


.0336 


± 





.3153 





.5235 


± 





.0731 


2 


,5102 


± 


0, 


,2791 


1 


.7313 


± 





,2625 





.7825 


± 





.1343 





9488 


± 





1568 


8.9 


4 


.2500 


± 


0. 


.1426 


1 


.2642 


± 





.0785 


2 


.9857 


± 





.1575 


2 


.7792 


± 





.2663 





.5656 


± 





.0578 


2 


,2136 


± 


0, 


,2340 


1 


,4708 


± 





,2005 





,6987 


± 





.0876 





,7721 


± 





1348 


9.0 


3 


7151 


± 





1449 


1 


.0853 


± 





.0615 


2 


6298 


± 





.1556 


2 


.4177 


± 





.2373 





.4832 


± 





.0559 


1 


9345 


± 


0, 


.1985 


1 


.2974 


± 


0, 


1458 





.6021 


± 





,0827 


0, 


6953 


± 





0795 


9.1 


3 


4459 


± 





.1257 


1 


.0891 


± 





.0444 


2 


.3568 


± 





1293 


2 


2779 


± 





.2120 





.5214 


± 





.0585 


1 


,7565 


± 





1743 


1 


1680 


± 


0, 


,1280 





.5677 


± 





.0704 





,6003 


± 





.0713 


9.2 


3 


.2127 


± 


0. 


1361 


1 


.0609 


± 





.0291 


2 


.1518 


± 





.1267 


2 


.1378 


± 





.1763 





.5271 


± 





.0433 


1 


.6106 


± 


0, 


1480 


1 


0749 


± 


0, 


0778 





.5337 


± 





.0439 


0, 


.5412 


± 





0465 


9.3 


2 


8792 


± 


0. 


.0998 


1 


.0113 


± 





.0308 


1 


8679 


± 





0993 


1 


9424 


± 





.1626 





.5276 


± 





.0440 


1 


.4148 


± 


0, 


1316 





9367 


± 


0, 


0924 





.4836 


± 





.0505 


0, 


,4531 


± 





0507 


9.4 


2 


.7471 


± 


0. 


0998 


1 


.0501 


± 





.0278 


1 


6970 


± 





0879 


1 


8702 


± 





.1530 





.5831 


± 





.0500 


1 


.2870 


± 


0, 


1125 





8770 


± 


0, 


0756 





.4670 


± 





.0434 


0, 


,4100 


± 





0393 


9.5 


2 


.6091 


± 


0. 


.0972 


1 


.0723 


± 





.0296 


1 


.5368 


± 





0783 


1 


.8150 


± 





.1525 





.6253 


± 





.0520 


1 


,1897 


± 


0, 


,1094 





7941 


± 





,0729 





,4470 


± 





.0362 





,3471 


± 





,0430 


9.6 


2 


6861 


± 





0956 


1 


.2200 


± 





.0332 


1 


4660 


± 





.0716 


1 


8920 


± 





.1553 





.7424 


± 





.0642 


1 


1496 


± 


0, 


0988 


0, 


.7940 


± 


0, 


.0747 





,4776 


± 





.0433 


0, 


.3164 


± 





0368 


9.7 


2 


.4858 


± 





.0875 


1 


.2074 


± 





.0339 


1 


.2784 


± 





0616 


1 


.7868 


± 





.1444 





.7688 


± 





.0662 


1 


,0180 


± 





.0856 


0, 


,6990 


± 





,0697 





,4386 


± 





.0419 


0, 


,2604 


± 





0315 


9.8 


2 


.1931 


± 





.0757 


1 


.1136 


± 





.0334 


1 


.0795 


± 





.0491 


1 


.6238 


± 





.1239 





.7524 


± 





.0579 


0, 


.8713 


± 





.0715 





,5693 


± 





,0564 





.3611 


± 





.0320 


0, 


,2082 


± 





0279 


9.9 


1 


.8351 


± 


0. 


0654 





.9485 


± 





.0274 


0. 


8866 


± 


0. 


0428 


1 


4036 


± 





.1038 





.6725 


± 





.0489 





.7311 


± 


0, 


.0593 





.4315 


± 


0, 


0450 





,2760 


± 





,0269 


0, 


1555 


± 





0205 


10.0 


1 


.5384 


± 


0. 


0557 





.8246 


± 





.0255 


0. 


.7138 


± 





0340 


1. 


.1979 


± 





0884 





.6080 


± 





.0452 





5899 


± 


0, 


.0465 


0, 


3405 


± 


0, 


0383 





,2165 


± 





.0241 


0, 


1239 


± 





.0160 


10.1 


1 


.2460 


± 





.0444 





.6707 


± 





.0210 


0. 


.5753 


± 





0265 





.9992 


± 





.0702 





.5149 


± 





.0364 





,4843 


± 


0, 


,0366 





2468 


± 





,0297 





.1558 


± 





.0186 





,0910 


± 





,0124 


10.2 





.9413 


± 





.0320 





.5112 


± 





.0160 


0. 


.4301 


± 





0182 





.7774 


± 





.0525 





.4108 


± 





.0274 





,3666 


± 


0, 


,0272 





1639 


± 





,0229 





.1004 


± 





.0133 





,0635 


± 





,0105 


10.3 





6428 


± 





0226 





.3591 


± 





.0123 





2837 


± 





.0120 





5448 


± 





0353 





.3000 


± 





.0190 





,2448 


± 


0, 


,0177 





.0980 


± 





,0141 





.0591 


± 





,0080 


0, 


0389 


± 





.0067 


10.4 





3961 


± 





.0137 





.2302 


± 





.0084 





.1659 


± 





0065 


0. 


.3442 


± 





0209 





,1995 


± 





,0121 





,1447 


± 





,0098 





0520 


± 





,0081 





,0307 


± 





,0046 


0, 


,0212 


± 





0039 


10.5 





.2227 


± 





.0075 





.1393 


± 





.0053 





.0834 


± 





.0030 


0. 


.1981 


± 





0104 





.1241 


± 





,0067 





0740 


± 





,0044 





,0246 


± 





,0035 





.0153 


± 





.0020 





,0094 


± 





0018 


10.6 





.1078 


± 





0028 





.0733 


± 





.0030 


0. 


0345 


± 





.0010 


0. 


0987 


± 





0042 





.0674 


± 





.0036 





0313 


± 





.0011 





0091 


± 





0018 





.0059 


± 





,0009 


0, 


0032 


± 





,0011 


10.7 





0475 


± 


0. 


0010 





.0345 


± 





.0012 


0. 


0130 


± 





.0011 





0447 


± 





0012 





.0328 


± 





.0013 





.0119 


± 





,0006 





0027 


± 





,0008 





.0017 


± 





.0003 


0, 


0011 


± 





,0006 


10.8 





.0183 


± 





.0011 





.0140 


± 





.0006 


0. 


.0043 


± 





.0010 





.0175 


± 





.0010 





.0135 


± 





,0006 





,0040 


± 





,0009 





0008 


± 





,0003 





,0005 


± 





.0001 





,0003 


± 





,0002 


10.9 





0060 


± 





.0006 





.0047 


± 





.0003 





0013 


± 





.0004 





.0059 


± 





.0005 





.0046 


± 





,0003 





0013 


± 





,0004 





.0001 


± 





,0001 





,0000 


± 





.0000 


0, 


.0000 


± 





,0001 


11.0 





0013 


± 





.0002 





.0011 


± 





.0002 





0002 


± 





.0001 


0. 


.0013 


± 





.0002 





.0011 


± 





,0002 





0002 


± 





,0001 





.0000 


± 





,0000 





,0000 


± 





,0000 


0, 


.0000 


± 





,0000 



Note. — Column (1): the median of the logarithm of the galaxy luminosity with bin width A log L — 0.05. Columns (2 - 4): the luminosity functions of all, red and blue galaxies for 'ALL' 
group members. Columns (5 - 7): the luminosity functions for all, red and blue galaxies for 'CENTRAL' group members. Columns (8 - 10): the luminosity functions of all, red and blue galaxies 
for 'SATELLITE' group members, respectively. Note that all the galaxy luminosity functions listed in this table are in units of 10 -2 /i 3 Mpc~ 3 d log L, where log is the 10 based logarithm. 
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TABLE 2 

The galaxy stellar mass functions $(M«) 



ALL CENTRAL SATELLITE 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 

logM* all red blue all red blue all red blue 



8.2 


6 


.7578 


± 


3 


5635 


5 


0225 


± 


3 


4088 


1 


7354 


± 


1 


3697 


5. 


4097 


± 


3 


6339 


3 


.6743 


± 


3 


1606 


1 


7354 


± 


1 


.3697 


1. 


3481 


± 


1 


4459 


1 


.3481 


± 


1. 


4459 





0000 


± 


o 


.0000 


8.3 


5 


.0065 


± 


1 


.1628 


2 


.6774 


± 


0. 


.9420 


2 


3291 


± 


o 


8917 


2 


4939 


± 





7368 


1 


.3292 


± 





6210 


1 


.1647 


± 


o 


.6085 


2 


.5126 


± 


o 


.7970 


1 


.3483 


± 


0. 


5956 


1 


.1644 


± 





5844 


8.4 


5 


2026 


± 





8331 


1 


.8689 


± 


0. 


5549 


3 


3337 


± 


o 


5824 


3 


7970 


± 





8574 





.9167 


± 





3865 


2 


8803 


± 


o 


.6688 


1 


4056 


± 


o 


.7744 





.9522 


± 





5715 





4534 


± 





.3364 


8.5 


3 


6943 


± 





6310 


1 


3602 


± 


0. 


3937 


2. 


3340 


± 


o 


.4160 


2 


8124 


± 





.5254 





.8612 


± 





2399 


1 


9512 


± 


o 


.4411 





.8819 


± 


o 


.5271 


o 


.4990 


± 





3365 





3829 


± 





.2879 


8.6 


3 


.2254 


± 





.4123 





.9867 


± 


0. 


2701 


2 


.2387 


± 


o 


3139 


2 


.4063 


± 





.4036 


o 


.5364 


± 


0. 


.1443 


1 


.8699 


± 





.3612 





8191 


± 


o 


.3452 


o 


4503 


± 


0. 


.2265 





.3687 


± 





.1858 


8.7 


3 


.1947 


± 





.3411 





8821 


± 


0. 


2339 


2 


.3126 


± 





.2237 


2 


.0724 


± 





.3095 


o 


.4822 


± 


0. 


.1251 


1 


5902 


± 


o 


.2715 


1 


.1223 


± 


o 


.3426 


o 


.3999 


± 





.1989 





.7224 


± 





.2100 


8.8 


2 


.7979 


± 





.2362 


o 


8262 


± 


0. 


1388 


1 


.9717 


± 


o 


2255 


2 


0071 


± 





.2667 


o 


.3576 


± 





.0824 


1 


.6495 


± 


o 


.2453 


o 


7908 


± 


o 


.1654 


o 


.4687 


± 





.1130 





3222 


± 





.0997 


8.9 


2 


7209 


± 


o 


.1921 


o 


.9936 


± 





1432 


1 


.7273 


± 





1742 


1 


.8136 


± 





2374 


o 


.4987 


± 


0. 


.0766 


1 


.3149 


± 


o 


.2136 


o 


9073 


± 


o 


.1880 


o 


.4949 


± 





.1142 


o 


.4124 


± 


o 


1215 


9.0 


2 


5366 


± 


o 


.1544 


o 


.7511 


± 





.1110 


1 


.7855 


± 





1573 


1 


.7138 


± 





1638 


o 


.3434 


± 


0. 


.0637 


1 


.3704 


± 


o 


.1687 





8228 


± 


o 


.1219 


o 


.4076 


± 


0. 


.0788 


o 


.4151 


± 





.0770 


9.1 


2 


8001 


± 


o 


.1436 


o 


.9423 


± 





1323 


1 


8578 


± 





.1762 


1 


8134 


± 





.2080 


o 


.3792 


± 





.0744 


1 


4342 


± 


o 


.1698 





.9867 


± 


o 


.1670 


o 


5630 


± 





.1391 





4236 


± 





.0625 


9.2 


2 


.5299 


± 





.1329 


o 


8579 


± 





1015 


1 


.6720 


± 





1546 


1 


.6551 


± 





.1690 


o 


.3732 


± 


0. 


.0545 


1 


2819 


± 





.1469 


o 


8748 


± 


o 


.1109 


o 


.4847 


± 


0. 


.0974 


o 


.3901 


± 





.0431 


9.3 


2 


.6046 


± 





.1105 


o 


8839 


± 


0. 


.0824 


1 


.7207 


± 


o 


1232 


1 


.7112 


± 





.1675 


o 


.4203 


± 


0. 


.0540 


1 


2908 


± 





.1388 





8935 


± 


o 


.1283 


o 


.4636 


± 





.0881 





.4299 


± 





.0612 


9.4 


2 


6018 


± 





.1118 


1. 


.0060 


± 





.0732 


1 


.5958 


± 





.1207 


1 


6445 


± 





.1602 





.4455 


± 





.0551 


1 


1990 


± 





.1283 





.9573 


± 





.1105 





.5605 


± 


0. 


.0818 





3968 


± 





.0480 


9.5 


2 


.2135 


± 





.1075 





.8682 


± 


0. 


.0492 


1 


.3453 


± 





.1048 


1 


4016 


± 





.1165 





.4033 


± 





.0414 


0. 


.9983 


± 





.0926 





.8118 


± 





.0559 





.4649 


± 





.0533 


0. 


3470 


± 





.0317 


9.6 


2. 


.1244 


± 





.0969 





.9046 


± 





.0347 


1 


.2197 


± 





.0840 


1 


.4653 


± 





1207 





.4732 


± 


0. 


.0443 





9921 


± 





.0951 





.6590 


± 





.0572 





.4314 


± 





.0407 


0. 


.2276 


± 





.0288 


9.7 


1 


.9695 


± 





.0934 





.9288 


± 





.0355 


1 


0407 


± 





.0762 


1 


.3411 


± 





.1129 





.5037 


± 





.0502 





8374 


± 





.0763 





6284 


± 





.0531 





.4251 


± 





.0434 


0. 


2033 


± 





.0194 


9.8 


1 


.7968 


± 





.0687 





.8633 


± 





.0252 


0. 


9335 


± 





.0614 


1 


2160 


± 





.0882 





.4676 


± 





.0395 





7483 


± 





.0616 





5808 


± 





.0451 





.3956 


± 


0. 


.0395 


0. 


.1852 


± 





.0144 


9.9 


1 


.7543 


± 





.0808 





.9275 


± 


0. 


.0410 





8268 


± 





.0538 


1 


1954 


± 





.0921 





.5415 


± 





.0491 





6539 


± 





.0545 





5589 


± 





.0340 





.3859 


± 


0. 


.0271 


0. 


.1729 


± 





.0134 


10.0 


1 


.6529 


± 





.0658 





.9104 


± 


0. 


.0335 





.7426 


± 





.0427 


1 


.1378 


± 





.0845 





.5414 


± 


0. 


.0472 





.5963 


± 





.0451 





.5152 


± 





.0338 





.3689 


± 


0. 


.0283 





.1462 


± 





.0115 


10.1 


1 


.6162 


± 





.0660 





.9446 


± 





0356 





.6716 


± 





.0398 


1 


.1427 


± 





.0837 





.5939 


± 





.0483 





5488 


± 





.0431 





4735 


± 





.0316 





.3507 


± 





.0264 





.1228 


± 





.0100 


10.2 


1 


.5700 


± 





.0631 





.9791 


± 





0388 





.5908 


± 





.0340 


1 


.1120 


± 





.0814 





.6402 


± 


0. 


.0531 





.4718 


± 





.0351 





.4580 


± 





.0310 





.3389 


± 





.0276 





.1190 


± 





.0076 


10.3 


1 


5025 


± 





.0564 





.9766 


± 


0. 


.0325 





.5259 


± 





.0298 


1 


.0918 


± 





.0766 





.6635 


± 


0. 


.0503 





.4284 


± 





.0320 





.4107 


± 





.0308 





.3132 


± 





.0271 





0975 


± 





.0066 


10.4 


1 


3108 


± 





.0480 





.8909 


± 


0. 


.0302 


0. 


4199 


± 





.0230 


0. 


9663 


± 





.0682 





.6177 


± 





.0479 





3486 


± 





.0252 





.3445 


± 





.0285 





.2732 


± 





.0247 


0. 


.0713 


± 





.0057 


10.5 


1 


0826 


± 





.0356 





.7608 


± 





0245 


0. 


3218 


± 





.0155 


0. 


8213 


± 





.0529 





.5522 


± 





.0384 





.2691 


± 





.0186 





.2613 


± 





.0237 





.2087 


± 


0. 


.0197 


0. 


0527 


± 





.0057 


10.6 





.8499 


± 





.0267 





.6080 


± 


0. 


.0188 





.2419 


± 





.0111 





.6604 


± 





.0418 





.4549 


± 


0. 


.0309 





2055 


± 





.0138 





.1896 


± 





.0193 





.1532 


± 





.0159 





.0364 


± 





.0046 


10.7 





.6440 


± 





.0193 





.4772 


± 


0. 


.0145 





.1667 


± 





.0074 


0. 


.5182 


± 





.0311 





.3753 


± 





.0242 


0. 


.1429 


± 





.0092 





.1257 


± 





.0151 





.1019 


± 


0. 


.0129 


0. 


0238 


± 





.0032 


10.8 





.4582 


± 





.0130 





.3556 


± 





0107 


0. 


1026 


± 





.0039 


0. 


.3811 


± 





.0211 





.2910 


± 





.0174 





.0901 


± 





.0051 





.0771 


± 





.0101 





.0646 


± 





.0087 





0125 


± 





.0020 


10.9 





3009 


± 





.0071 





.2436 


± 





.0069 





.0573 


± 





.0014 





.2572 


± 





.0123 





.2064 


± 


0. 


.0114 





.0508 


± 





.0019 





0437 


± 





.0067 





.0373 


± 





.0059 





.0064 


± 





.0012 


11.0 





1846 


± 





.0047 





.1573 


± 


0. 


.0048 





.0273 


± 





.0011 





.1627 


± 





.0073 





.1379 


± 


0. 


.0071 





.0248 


± 





.0008 





0219 


± 





.0035 





.0194 


± 





.0031 





.0025 


± 





.0006 


11.1 





1066 


± 





.0022 





.0931 


± 





.0027 


0. 


0135 


± 





.0012 


0. 


0964 


± 





.0037 





.0839 


± 





.0039 





0125 


± 





.0008 





0102 


± 





.0020 





.0092 


± 





.0016 


0. 


0010 


± 





.0005 


11.2 





0575 


± 





.0011 





.0520 


± 





.0017 


0. 


0055 


± 





.0009 


0. 


0531 


± 





.0018 





.0480 


± 





.0023 


0. 


.0052 


± 





.0007 





.0044 


± 





.0010 





.0010 


± 


0. 


.0008 


0. 


0003 


± 





.0002 


11.3 





0294 


± 





.0006 





.0271 


± 


0. 


.0009 





.0023 


± 





.0004 


0. 


.0281 


± 





.0009 





.0258 


± 


0. 


.0011 


0. 


0022 


± 





.0004 





.0014 


± 





.0004 





.0013 


± 


0. 


.0003 


0. 


.0001 


± 





.0001 


11.4 





.0131 


± 





.0004 





.0122 


± 


0. 


.0005 





.0010 


± 





.0002 


0. 


.0127 


± 





.0005 





.0117 


± 





.0005 





.0009 


± 





.0002 





0004 


± 





.0001 





.0004 


± 


0. 


.0001 


0. 


0000 


± 





.0000 


11.5 





.0047 


± 





.0003 





.0044 


± 





0003 


0. 


0003 


± 





.0001 


0. 


0047 


± 





.0003 





.0044 


± 





.0003 





.0003 


± 





.0001 





0000 


± 





.0000 





.0000 


± 





.0000 





.0000 


± 





.0000 


11.6 





.0012 


± 





.0001 





.0011 


± 





.0001 





.0001 


± 





.0000 


0. 


.0012 


± 





.0001 





.0011 


± 


0. 


.0001 





.0001 


± 





.0000 


0. 


0000 


± 





.0000 





.0000 


± 





.0000 





.0000 


± 





.0000 



Note. — Column (1): the median of the logarithm of the galaxy stellar mass with bin width A log Af* — 0.05. Column (2 - 4): the stellar mass functions of all, red and blue for 'ALL' group 
members. Column (5 - 7): the stellar mass functions of all, red and blue for 'CENTRAL' group members. Column (8 - 10): the stellar mass functions of all, red and blue for 'SATELLITE' group 
members. Note that all the galaxy stellar mass functions listed in this table are in units of 10 -2 ft, 3 Mpc -3 d log M* , where log is the 10 based logarithm. 
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Fig. 2. — The galaxy luminosity functions (left panels) and stellar mass functions (right panels). The upper, middle and lower panels 
show results obtained from all, central and satellite group members, respectively. In each panel, the open circles, squares and triangles 
with error-bars are the luminosity or stellar mass functions of all, red and blue galaxies, respectively, where the errors are obtained from 
the 200 bootstrap samples. Note that the results for red and blue galaxies are scaled down by a factor of 10, for clarity. The solid lines in 
each panel are the best fitting Schechter functions. For comparison, the dot-dashed line in the upper-left panel is the best fit luminosity 
function obtained by Blanton et al. (2003b), while that in the upper-right panel corresponds to the stellar mass function obtained by Bell 
et al. (2003) from the SDSS Early Data Release. 



luminosity L* (or stellar mass M*). Using the least \ 2 
fitting to the measured luminosity and stellar mass func- 
tions shown in Fig. [2j we find the best fit values of </>*, 
a, L* (or M*), which are listed in Tabled The best fit 
results are shown as the solid lines in Fig.O 

Comparing the data with the best fit, one sees that the 
Schechter form describes the data remarkably well over 
the luminosity range logL > 9.0. At the fainter end, 
however, the data reveals an upturn, which is almost en- 
tirely due to the red centrals. The stellar mass functions 
also show a steepening at the low mass end, which again 
owes mainly to the population of red centrals. 

The galaxy luminosity and stellar mass functions have 
been estimated before by various authors. As an illustra- 
tion, we show the best fit of the luminosity function ob- 
tained by Blanton et al. (2003b) from the SDSS DR2 as 
the dot-dashed line in the upper-left panel of Fig. [2l As 
expected, our measurements are in excellent agreement 
with theirs. For the stellar mass function, we show as the 
dot-dashed line the result of Bell et al. (2003) obtained 



from the SDSS Early Data Release (EDR; Stoughton et 
al. 2002). The agreement with our measurement is re- 
markably good for logM* > 9.5. For lower masses, how- 
ever, the mass function obtained by Bell et al. (2003) is 
significantly higher. The discrepancy most likely results 
from the different data samples (DR4 vs. EDR) used in 
the two analyses. For low-mass galaxies, the cosmic vari- 
ance is significant, especially in EDR because only a few 
hundred galaxies in a small volume were used to measure 
the stellar mass function. In addition, we have taken into 
account the redshift completeness of galaxies using the 
completeness masks provided by the NYU team, and the 
stellar mass limit is treated more carefully in our analysis 
using Eq. [Jj The behavior of the stellar mass function in 
the low-mass end has also been investigated by Baldry et 
al. (2004, 2008) and Panter et al. (2007). Unfortunately, 
the situation is still unclear, partly because of the limited 
sample volume, and partly because of the uncertainties 
in the luminosity-mass conversion. 
To see where galaxies of different luminosities and 
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Fig. 3. — The host halo mass distributions for central (red lines) and satellite (blue lines) galaxies. Shown in the upper panels are results 
for galaxies in different luminosity bins as indicated. Shown in the lower panels are results for galaxies in different stellar mass bins as 
indicated. The dotted and solid lines represent the results for halo masses, Mr, and Ms, respectively. 



TABLE 3 

The best fit parameters for the galaxy 
luminosity functions and stellar mass 
functions 



Member type 


color 




a 


logL* 


(1) 


(2) 


(3) 


(4) 


(5) 


ALL 


all 


0.03167 


-1.117 


10.095 




red 


0.01810 


-0.846 


10.080 




blue 


0.01890 


-1.154 


10.026 


CENTRAL 


all 


0.02370 


-1.069 


10.114 




red 


0.01223 


-0.755 


10.106 




blue 


0.01523 


-1.123 


10.033 


SATELLITE 


all 


0.01051 


-1.134 


9.947 




red 


0.00682 


-1.018 


9.923 




blue 


0.00410 


-1.235 


9.951 


Member type 


color 




a 


log M* 


ALL 


all 


0.01546 


-1.164 


10.717 




red 


0.01285 


-0.916 


10.691 




blue 


0.00836 


-1.233 


10.508 


CENTRAL 


all 


0.01084 


-1.143 


10.758 




red 


0.00903 


-0.803 


10.717 




blue 


0.00668 


-1.218 


10.522 


SATELLITE 


all 


0.00692 


-1.078 


10.483 




red 


0.00589 


-0.932 


10.457 




blue 


0.00165 


-1.290 


10.450 



Note. — Column (1): the member type. Col- 
umn (2): the color of galaxies. Column (3-5): the 
best fit parameters for the luminosity functions (up- 
per part) and stellar mass functions (lower part). Note 
that <fi* listed in column 3 are presented in terms of 
7i 3 Mpc~ 3 dlogL (or /i 3 Mpc~ 3 d log M, ), where log is 
the 10 based logarithm. 

stellar masses are hosted, we plot in Fig. [3] the host 
halo mass distribution for central (red lines) and satel- 
lite (blue lines) galaxies. Results are shown for galaxies 
in different luminosity (upper panels) and stellar mass 
(lower panels) bins, as indicated 9 . Almost all bright (or 
massive) galaxies are centrals in massive halos. About 

9 Note that in the group catalog, halo masses are provided only 



3/4 of faint (low-mass) galaxies are centrals in small ha- 
los with a very narrow mass distribution, and the rest 1/4 
are satellites in halos that cover a wide range in mass. 
The results are shown separately for halo masses, Ml 
and Ms- Although the results for satellite galaxies and 
the mean halo mass for central galaxies are similar for 
the two halo mass estimates, the widths of the halo mass 
distribution for central galaxies are quite different. This 
is caused by the fact that these is some spurious correla- 
tion between L19.5 (Medlar) and the luminosity (stellar 
mass) of the central galaxy, especially for low-mass ha- 
los where the luminosity (stellar) content is dominated 
by the central galaxy. We have measured the halo mass 
distribution from both samples II and III, and we do not 
find any significant difference between the results. There- 
fore only the results for sample II are plotted in Fig|3l 
The general behavior of the halo mass distribution of the 
central and satellite galaxies is consistent with that pre- 
dicted by the CLF and HOD models (e.g. Yang et al. 
2003; Zheng et al. 2005). 

4. THE CONDITIONAL STELLAR MASS 
FUNCTION 

In paper II, we have measured the conditional luminos- 
ity function (CLF) of galaxies in halos (as represented 
by galaxy groups), $(L|M,j). Here we first obtain the 
conditional stellar mass function (CSMF) of galaxies in 
dark halos. The CSMF, $(M*|M/,), which describes 
the average number of galaxies as a function of galaxy 
stellar mass in a dark matter halo of a given mass, is 
more straightforwardly related to theoretical predictions 
of galaxy formation models than the CLF, because the 
conversion from stellar mass to luminosity in theoretical 
models requires detailed modeling of the stellar popula- 
tion and dust extinction. The CSMF can be estimated 

for groups whose central galaxies are brighter than 01 M r — 5 log h = 
— 19.5. For group s w ith a fainter central galaxy we us e the mean 
mass-to-light (Eq[l9j and halo-to-stellar mass (Eg 1200 ratios ob- 
tained in Section [6] to estimate their halo masses. 




O 
u 

DO 



DO 
O 



o 



- 1 — i — i — i — i — i — i — i — i — i — i — i — i — r 

12.9<log M h ^13.2 




J I I I L 



9 10 11 12 
log[M„/(h-»M e )] 



o 
s- 



00 
O 



IS 

""So 
o 




9 10 11 
log[M„/(h- 8 M @ )] 



Fig. 4. — The conditional stellar mass functions (CSMFs) of galaxies in groups of different mass bins. Symbols correspond to the 
CSMFs obtained using Ms as halo mass (estimated according to the ranking of the characteristic group stellar masses), with solid and 
open circles indicating the contributions from central and satellite galaxies, respectively. The error-bars reflect the 1-a scat ter obt aine d 
from 200 bootstrap samples for Ml and M§, respectively. The solid lines indicate the best-fit parameterizations (equations 1131 to 1161 ). 
For comparison, we also show, with dashed lines, the CSMFs obtained using Ml as halo mass (estimated according to the ranking of the 
characteristic group luminosity). 



by directly counting the number of galaxies in groups. 
However, because of the completeness limits discussed 
in Section [321 we only use galaxies and groups that are 
complete according to Eqs. ([7]) and © to estimate the 
CSMF, $(M*\M h ), at a given M*. In Fig.[4]we show the 
resulting CSMFs for groups of different masses. The con- 
tributions of central and satellite galaxies are plotted sep- 
arately. For comparison, results obtained using Ms and 
Ml are shown as symbols and dashed lines, respectively. 
The error-bars shown in each panel correspond to 1-a 
scatter obtained from 200 bootstrap samples of our group 
catalogue. In general, these two halo masses give con- 
sistent results, except that the Ms-based CSMF of the 
central galaxies in low mass halos is more peaked than 
the Mi-based CSMF (see the lower right-hand panel). 
The general behavior of the CSMF is similar to that of 
the CLF presented in Paper II. The general behavior of 



the CSMF obtained here is also qualitatively similar to 
the prediction of semi-analytical models (e.g. Zheng et 
al. 2005): the CSMF for small halos has a strong peak 
at the bright end due to central galaxies. Quantitatively, 
however, semi-analytical models in general over-predict 
the number of satellite galaxies (Liu et al. in prepara- 
tion). 

In Fig. [5] we show the CSMFs separately for red 
(dashed lines with errorbars) and blue (dotted lines) 
galaxies. Clearly there are more red galaxies than blue 
galaxies (both centrals and satellites) in massive halos. 
In the lowest mass bin probed here (12.0 < logM^, < 
12.3), however, there are roughly equal numbers of red 
and blue galaxies. The fraction of red galaxies as a func- 
tion of halo mass found here is very similar to that ob- 
tained by Zandivarez et al. (2006) based on the condi- 
tional luminosity function of galaxies derived from an in- 
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Fig. 5. — Similar to Fig. [4] but here we show the CSMFs for red (dashed lines) and blue (dotted lines) galaxies. In both cases the central 
and satellite components of the CSMFs are indicated separately. The error-bars, again obtained using 200 bootstrap samples, are shown 
for the red and blue galaxies, separately. The solid lines indicate the best-fit parameterizations (equations 1131 to 1161 ) for red galaxies. 
Results shown are for halo masses Mg only. 



dependent group catalog. Note that the overall shapes of 
the CSMFs for red and blue galaxies are remarkably sim- 
ilar. Interestingly, such behavior is predicted by Skibba 
et al. (2008) who used the color-marked correlation func- 
tion to constrain the distribution of galaxies according to 
their colors. 

We model the CSMF using the sum of the CSMFs 
of central and satellite galaxies (see Yang et al. 2003; 
Cooray 2005; White et al. 2007; Zheng et al. 2007; Yang 
et al. 2008; Cacciato et al. 2008): 

$(M4M h ) = $ cea (M*\M h ) + $ sat (M„|M h ) . (13) 

Following Paper II, we adopt a lognormal model for the 
CSMF of central galaxies: 



$ cen (M*|i\4) 



.4 



27TCT, 



-exp 



(logM* - logM, 



2(7? 



(14) 

where A is the number of central galaxies per halo. Thus, 
A = 1 for all galaxies, A = f ro d(Mh) for red galaxies, and 



A = 1 — f rc d(M h ) for blue galaxies. Here f re d(M h ) is the 
red fraction of central galaxies in halos of mass Mh ■ Note 
that logM^c is, by definition, the expectation value for 
the (10-based) logarithm of the stellar mass of the central 
galaxy: 



logM* 



$ CC „(M, |M ft ) logM.dlogM, , (15) 



and that a c = <r(log JVf*). For the contribution from the 
satellite galaxies we adopt a modified Schechter function: 



$ sat (M*|M h ) 



\M* 



exp 



\M* 



(16) 

Note that this function decreases faster at the bright end 
than a Schechter function and gives a better description 
of the data. The above parameterization has a total of 
five free parameters: M* )C , a c , </>*, a* and M* iS . Wc 
find that logM*, c = logM* iS + 0.25 to good approxima- 
tion, which is what we adopt throughout. Consequently, 
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Fig. 6. — The best fit parameters ((f)* upper row, a* second row, M* :C third row, and <x c bottom row) to the CSMFs shown in Figs.[3]and 
[5] as functions of halo mass. Panels on the left, in the middle, and on the right show results for all, red, and blue galaxies, respectively. 
Since we have two different halo mass estimators (Mg and Ml) and two main group samples (II and III), we have obtained CLFs for four 
different combinations of sample and group mass estimator. The results for all four combinations are shown using different symbols and 
line-styles, as indicated. The error-bars in the first two and last rows indicate the l-tr variances obtained from our 200 bootstrap samples. 
In the third row of panels, however, the error-bars correspond to the log-normal scatter, <r c , shown in the bottom row of panels. For clarity 
the error-bars are only shown for the 'Mg-sample IF case, but they are very similar for the other four cases. 



the number of free parameters is reduced to four. Note, 
however, Hansen et al. (2007) found that the ratio be- 
tween the mean stellar mass (luminosity) of the central 
galaxy and the characteristic stellar mass (luminosity) of 
the satellite galaxies depends on halo mass especially for 
massive halos. The difference may arise from the fact 
that they are using galaxy groups constructed from the 
SDSS photometric redshift catalogue where the member- 
ships of the groups are not so well constrained. 

We fit the above model to all our CSMFs. Results are 
shown in Fig [51 separately for all (left panels), red (mid- 
dle panels) and blue (right panels) galaxies. Here com- 
parison is made between samples II and III. Note that 
Sample II does not include any galaxies missed due to 
fiber collisions, while Sample III includes all such galax- 
ies by assigning each of them the redshift of its nearest 
neighbor. Given that we also have two kinds of halo 
masses, Ml and Ms, there are a total of four different 
combinations for which we have determined the CSMFs. 
Each panel of Fig [5] shows the results for all four cases. 
As an illustration of how well the model fits the data, 



the solid lines in Figs. [Hand [5] show the corresponding 
best-fit models. 

The upper row of Fig. [6] shows the best fit normal- 
ization of the CSMF for satellite galaxies, which de- 
scribes the average number of satellite galaxies with stel- 
lar mass ~ M# tS in a group of a given halo mass. As 
expected, Sample III gives a higher </>*, especially for 
low-mass groups, but only marginally so. Comparing 
4>* for red (upper middle panel) and blue (upper right 
panel) galaxies, one sees that the fraction of red satel- 
lites increases with halo mass. The second row shows the 
faint end slopes of the CSMFs, a*. In massive halos with 
Mh > 1O 13 /i _1 M0, a* s decreases (i.e., becomes more neg- 
ative) with increasing halo mass, both for red and blue 
galaxies. In halos with Mh < 10 13 /i _1 M0, however, a* 
is roughly constant at ~ —1.2. The third row shows that 
logM* jC increases with halo mass, for both red and blue 
centrals. Note that for a given value of logM* jC , the halo 
mass for blue galaxies based on Ml is larger than that 
based on Mg, especially in small halos. The reason is 
that for a given stellar mass, bluer galaxies are brighter, 
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TABLE 4 

The best fit parameters of the CSMFs for all, red and blue galaxies 



Galaxy type 


log[Af h ] 


log([M h ]) 






logM», c 


Cc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




14.40, 15.00) 


14.58 


24.91 ±2.51 


-1.59 ±0.05 


11.364 ±0.018 


0.159 ±0.017 




14.10, 14.40) 


14.23 


16.49 ±0.79 


-1.40 ±0.03 


11.277 ±0.013 


0.164 ±0.011 




13.80, 14.10) 


13.94 


9.72 ± 0.63 


-1.34 ±0.07 


11.209 ±0.020 


0.182 ±0.011 




13.50, 13.80) 


13.64 


6.15 ±0.77 


-1.22 ±0.07 


11.122 ±0.020 


0.180 ±0.007 


ALL 


13.20, 13.50) 


13.35 


3.61 ±0.52 


-1.16 ±0.08 


11.026 ±0.024 


0.179 ±0.018 




12.90, 13.20) 


13.05 


1.96 ±0.47 


-1.15 ±0.14 


10.926 ± 0.032 


0.168 ±0.025 




12.60, 12.90) 


12.75 


1.10 ±0.22 


-1.19 ±0.09 


10.803 ± 0.034 


0.142 ±0.035 




12.30, 12.60) 


12.45 


0.61 ±0.10 


-1.19 ±0.06 


10.660 ± 0.033 


0.140 ±0.046 




12.00, 12.30) 


12.16 


0.31 ±0.06 


-1.24 ±0.09 


10.485 ± 0.022 


0.145 ±0.058 



RED 



14.40 
14.10 
13.80 
13.50 
13.20 
12.90 
12.60 
12.30 
12.00 



15.00) 
14.40) 
14.10) 
13.80) 
13.50) 
13.20) 
12.90) 
12.60) 
12.30) 



14.58 
14.23 
13.94 
13.64 
13.35 
13.05 
12.75 
12.45 
12.16 



22.64 ±2.66 
14.40 ± 0.74 
8.30 ±0.33 
4.86 ±0.73 
2.71 ±0.39 
1.37 ±0.24 
0.71 ±0.15 
0.36 ±0.09 
0.18 ±0.04 



-1.54 ±0.05 
-1.39 ±0.03 
-1.30 ±0.07 
-1.18 ±0.06 
-1.10 ±0.06 
-1.11 ±0.07 
-1.16 ±0.10 
-1.14 ±0.07 
-1.17 ±0.12 



11.365 
11.278 
11.216 
11.131 
11.047 
10.958 
10.844 
10.713 
10.539 



± 0.020 
±0.013 
±0.017 
±0.015 
±0.014 
±0.012 
± 0.009 
±0.027 
± 0.042 



0.162 
0.165 
0.179 
0.179 
0.175 
0.161 
0.127 
0.112 
0.100 



±0.018 
±0.012 
± 0.008 
± 0.008 
±0.015 
± 0.020 
± 0.027 
±0.018 
± 0.007 



BLUE 



14.40 
14.10 
13.80 
13.50 
13.20 
12.90 
12.60 
12.30 
12.00 



15.00) 
14.40) 
14.10) 
13.80) 
13.50) 
13.20) 
12.90) 
12.60) 
12.30) 



14.58 
14.23 
13.94 
13.64 
13.35 
13.05 
12.75 
12.45 
12.16 



3.23 ±5.24 
2.11 ± 1.06 
1.67 ±0.27 
1.41 ±0.26 
0.99 ±0.26 
0.68 ±0.14 
0.43 ±0.05 
0.26 ±0.04 
0.15 ±0.03 



-1.71 ±0.31 
-1.48 ±0.18 
-1.41 ±0.09 
-1.26 ±0.08 
-1.23 ±0.12 
-1.12 ±0.10 
-1.13 ±0.07 
-1.21 ±0.10 
-1.25 ±0.14 



11.305 
11.276 
11.152 
11.082 
10.978 
10.872 
10.746 
10.601 
10.414 



±0.134 
± 0.073 
± 0.035 
± 0.040 
± 0.053 
± 0.067 
± 0.069 
± 0.085 
± 0.097 



0.183 
0.174 
0.197 
0.183 
0.176 
0.174 
0.135 
0.127 
0.122 



± 0.049 
± 0.044 
±0.013 
±0.012 
±0.017 
±0.010 
± 0.020 
± 0.027 
± 0.030 



Note. — Column (1): galaxy type. Column (2): halo mass range. Column (3): average of the logarithm of 
the halo mass. Column (4)-(7): average of the best fit free parameters to the four measurements of the CSMFs. 
as shown in Fig. \§\ The errors indicate the scatter among these four measurements or the scatter obtained from 
the 200 bootstrap samples, whichever is larger. 



hence Ml is higher. This effect can impact (e.g. slightly 
enhance) our results on the color dependence of group 
clustering (e.g., Berlind et al. 2006; Wang et al. 2008a). 
These best fit parameters for satellite galaxies of differ- 
ent colors are different from that obtained from the CLF 
measurements (e.g. Zandivarez et al. 2006; Hansen et al. 
2007; Yang et al. 2008). In particular the low-mass end 
slope, a* , obtained here for red satellite galaxies in small 
halos is significantly steeper. The only significant differ- 
ence between red and blue galaxies is in </>*, consistent 
with the results of Skibba (2008). Finally, the last row of 
Fig [6] shows the width of the log-normal CSMF of central 
galaxies. For the combined sample of red and blue galax- 
ies we find an average value of a c — o^logM*^) ~ 0.17 
in halos with M h > 10 13 /j _1 M Q . This width is slightly 
larger than that in the luminosity distribution obtained 
in Paper II. A roughly constant dispersion in the lognor- 
mal distribution of luminosity (or stellar mass) for central 
galaxies has already been predicted by Yang et al. (2003), 
Cooray (2006) and Cacciato et al. (2008) from the clus- 
tering and abundances of galaxies. Although Zheng et 
al. (2007) found, based on a HOD model applied to the 
SDSS and DEEP2, that the log-normal width increases 
from ~ 0.13 for massive halos with Mh ~ 10 13 5 /t _1 M 
to ~ 0.3 for low mass halos with M h ~ 10 115 hr l M Q , 
the difference between their results and ours is only at 
l-<7 level. Again, since our halo masses are based on the 
ranking of either £19.5 or Msteiiar, we have effectively as- 
sumed a one-to-one relation between halo mass and these 



mass indicators. This can give rise to spurious correla- 
tion between L19.5 (or A/ ste iiar) with the stellar mass of 
the central galaxies, especially for low-mass halos where 
the stellar content is dominated by the central galaxy. 
Therefore the values of er c obtained, especially for halos 
with Mh < 10 13 h~ 1 M Q , may be underestimated and 
should be considered as lower limits. 

For reference, Table [4] lists the average values of the 
CSMF fitting parameters obtained from the combina- 
tions of Samples II and III and group masses Ml and 
Ms- The error-bars indicate the scatter among these 
four samples or the scatter obtained from 200 bootstrap 
samples, whichever is larger (generally the former). 

5. THE GROUP LUMINOSITY AND 
STELLAR-MASS FUNCTIONS 

In this section, we present our results on the group lu- 
minosity and stellar-mass functions. These two functions 
depend the total luminosity and stellar mass in galaxy 
groups, and are arguably better suited for comparison 
with model predictions, as the details about how the 
total luminosity and total stellar mass are partitioned 
into individual member galaxies is not important here. 
We measure the group luminosity and stellar-mass func- 
tions for both samples II and III. As mentioned above, 
these two samples represent two extremes, as far as fiber- 
collision effects are concerned. The results presented in 
the following are based on the average of samples II and 
III (evenly weighted), and the error-bars are obtained 
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Fig. 7. — The group luminosity functions and stellar mass functions. The upper-left panel: solid and dotted histograms are the group 
luminosity functions for -L19.5 and I/total > respectively. The open circles with error-bars are the sum of the two contributions. The dashed 
line is the luminosity function for central galaxies shown in Fig. [2] The upper-right panel: solid and dotted histograms are the group stellar 
mass functions for iW^ellar an d Af tota i, respectively. The open circles with error-bars are the sum of the two contributions. The dashed 
line is the stellar mass function for central galaxies shown in Fig. [2] The lower-left panel: the open circles with error-bars are the same as 
those in the upper-left panel and the long dashed line shows the best fitting results. The solid histogram is the group luminosity function 
for 1/18.5- The lower-right panel: the open circles with error-bars are the same as those in the upper-right panel and the long dashed line 
shows the best fitting results. The solid histogram is the group stellar mass function for Mis. 5. See text for the definitions of L19.5, -Ltotab 
Lis. 5, Af stc n ar , M tota i and Mis. 5. For comparison, in the lower two panels, we show as dot-dashed lines the scaled halo mass function 
(Warren et al. 2006) for the same ACDM 'concordance' cosmology. 



from the difference between the two samples, or from 200 
bootstrap samples for samples II and III, whichever (typ- 
ically the difference between the two samples) is larger. 

As discussed in Section l2~2"l the group catalogue suffers 
from incompleteness, in that groups within certain lumi- 
nosity (L19.5), stellar mass (M stc iiar) and halo mass bins 
are only complete to certain redshifts. Using the same 
procedure as described in Section 12.21 we estimate, for 
each group luminosity Lig.5 (or stellar mass M st ciiar), the 
maximum redshift, z max , that satisfies Eq. JSJ. The re- 
sulting n(z max ) as a function of group luminosity £19.5 or 
stellar mass M st eiiar gives, the group luminosity function 
^(Lig.s) or the group stellar mass function $(M s teiiar)) 
respectively. The results are shown as the solid his- 
tograms in the upper panels of Fig. [7J 



Note that the group characteristic luminosity Lig.5 
and stellar mass M st ciiar are defined to be the total 
luminosity and stellar mass of member galaxies with 
01 M r — 5 log ft. < —19.5. There are certainly groups in 
which all member galaxies have 01 M r — 5 log h > —19.5, 
so that Lig.5 = and Miliar = 0. For these groups 
we measure their total luminosity itotai and total stel- 
lar mass Mtotai based on all member galaxies that are 
observed. We measure the luminosity function <E>(Ltotai) 
and stellar mass function $(M to tai) for these groups us- 
ing the same method as for <J>(Lig. 5 ) and ^(Mstciiar), 
i.e., by measuring the n(z m&x ). The results are shown 
in the upper panels of Fig [7] as the dotted histograms. 
For comparison, we also show as the dashed lines the 
galaxy luminosity function (stellar mass function) for 
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TABLE 5 

The group luminosity functions and stellar mass functions 



logL(logM«) $(£total) HLl9.5) ^(LlS.5) $(jjtotal) gWrteUar) ^(^18.5) 
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0. 
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.0000 
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Note. — Column (1): the median of the logarithm of the group luminosity or stellar mass with bin width A log L — 0.05 (or 
A log M* — 0.05 ). Column (2): the average group luminosity function for L tota i, where L tota i is the total luminosity of all group 
members in which the central galaxy has luminosity 01 M r — 5 log h > —19.5. Column (3): the average group luminosity function for 
L19.5, where Lig.5 is the luminosity of all group members with *M r — 5 log h < —19.5. Column (4) the average group luminosity 
function for Lis. 5, where Lis. 5 is the total luminosity of all group members with 01 M r — 5 log h < —18.5. Columns (5)-(7): Similar 
to Columns (2)-(4) but for the average stellar mass functions. In this table, the averages arc obtained from the two measurements of 
the luminosity functions or stellar mass functions from samples II and III, respectively. The error indicates the scatter among these 
two measurements or the scatter obtained from 200 bootstrap samples, whichever is larger. Note that all the group luminosity (stellar 
mass) functions listed in this table are calculated in units of 10 -2 /i 3 Mpc — 3 d log L (or 10 — 2 h 3 Mpc — 3 d log A-f* ), where log is the 10 based 
logarithm. 



central galaxies that we obtained in Section[3] The agree- 
ment between the group luminosity function $ (Ltotai) 
(stellar mass function ^(Mtotai)) and the central galaxy 
luminosity function (stellar mass function) at the faint 
(low-mass) end is excellent. This simply reflects that 
the total luminosity Ltotai and total stellar mass M to tai 
of small groups are dominated by their central galaxies 
(e.g. Lin & Mohr 2004; Hansen et al. 2007; Yang et al. 
2008b), and implies that missing faint satellites in low 
mass groups has no significant impact on our estimates 



of the group luminosity and stellar-mass functions, unless 
there is a sharp upturn in the conditional luminosity (or 
stellar mass) function of satellite galaxies, as indicated in 
some observations (e.g. Popesso et al. 2006). Thus the 
overall group luminosity function (stellar mass function) 
from small to large groups can be obtained by sum up 
$(Li 9 . 5 ) and $(L t otai), WMrteilar) and $(M tot ai))- The 
results are shown as the open circles in the upper panels. 

As a further test of the reliability of our results at 
the low-luminosity and low-mass end, we use another 
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characteristic group luminosity Lis. 5 and another char- 
acteristic stellar mass Mis. 5 to represent the luminosity 
and stellar mass of groups. Here Lis. 5 (M18.5) is de- 
fined to be the total luminosity (stellar mass) of all group 
members with 01 M r — 5 log ft, < —18.5. Because of the 
survey apparent magnitude limit, a galaxy fainter than 
01 M r — 5 log ft = —18.5 is not observed if its redshift is 
larger than ~ 0.06. To avoid such incompleteness, Lis. 5 
and Mis. 5 are only measured for groups with redshift 
z < 0.06. Our method to estimate ^(Lis.s) and ^(Mis.s) 
is the same as that used for $(Lig. 5 ) and $(M stc iiar), ex- 
cept that z m ax starts from 0.06 instead of z = 0.2. The 
results for $(Li 8 .5) and $(Mis.s) are shown in the lower 
two panels of Fig. [7] as the solid histograms. For compar- 
ison, we also show L19.5 and M s t e iiar in the corresponding 
panels as open circles. Clearly, there is a good agreement 
between Q(L\%.§) and <j>(Lig. 5 ), and between $(Mi8.s) 
and $ (Miliar) for small groups, suggesting again that 
our results for low-mass groups are reliable. For massive 
groups, there is some difference between the two esti- 
mates. Unfortunately, the number of such groups is quite 
small in the volume corresponding to z = 0.06, and the 
results for $(Lis.s) and $(Mi8.s) are quite uncertain at 
the massive end. For reference, the data for the group 
luminosity functions and stellar mass functions shown in 
Fig. [7] are listed in Table [5] 

Since both the luminosity function and stellar mass 
function for groups show double power-law behavior, we 
use the following forms to fit the data: 



$(L 



19.5; 



$0 



(Ll9.5/Lo) C 



(a* + (Lig.s/Lo) 4 )' 3 



and 



$(M 



stellar J 



$0- 



(M stcllar /M )" 
x + (M stcUar /M ) 4 )' 3 



(17) 



(18) 



The best fit results are shown in the lower two panels 
of Fig. [7j as the long-dashed lines. The corresponding 
best fit parameters are [$o, logLo, xo, a, /?] =[0.00580, 
10.30, 0.1786, -0.2226, 0.4236] and [$ , log M , x , a, 0\ 
= [0.00731, 10.67, 0.7243, -0.2229, 0.3874 ], for the lu- 
minosity and stellar-mass functions, respectively. Note 
that all the group luminosity functions and stellar mass 
functions are calculated in terms of ft 3 Mpc _3 d log L and 
ft 3 Mpc~ 3 dlogM i( ,), where log is the 10 based logarithm. 

It is interesting to compare the group luminosity and 
stellar-mass functions with the halo mass function pre- 
dicted by the current CDM model of structure forma- 
tion. In the lower-left panel of Fig we plot, as 
the dot-dashed line, the halo mass function which is 
scaled with a (minimum) constant mass-to-light ratio, 
M/j/Lig 5 = 47/iMq/L q , so that the halo mass function 
touches the group luminosity function at log L19.5 ~ 10.2. 
Here the halo mass function is estimated using the War- 
ren et al. (2006) model and adopting the ACDM 'con- 
cordance' cosmology with parameters listed at the end 
of Section 1. Clearly, the scaled halo mass function 
has a very different shape from the observed group lu- 
minosity function, indicating that the mass-to-light ra- 
tio must depend strongly on halo mass. The mass-to- 
light ratio reaches its minimum in groups (halos) with 
logLig.5 ~ 10.2, and is larger for both smaller and big- 
ger halos. In the lower-right panel of Fig. [7] we com- 
pare the group stellar mass function (open circles) with 



the halo mass function scaled with a constant halo- to- 
stellar mass ratio (dot-dashed line), Mh/M ste uar = 25ft. 
Here the halo-to-stellar mass ratio reaches its minimum 
in groups (halos) with logM st oiiar ~ 10.6, corresponding 
to halos with masses Mh ~ 11.9 (almost exactly the same 
halo mass that the group mass-to-light ratio reaches its 
minimum), suggesting that star formation efficiency is 
the highest in such halos. These results are in excel- 
lent agreement with those obtained from the clustering 
properties of galaxies using the CLF and HOD formalism 
(Yang et al. 2003; Tinker et al 2005; van den Bosch et al. 
2007; Cacciato et al. 2008), and at the relative massive 
end with those obtained from weak lensing observations 
(e.g. Sheldon et al. 2007). 

6. THE CENTRAL GALAXIES IN SMALL HALOS 

In this section we discuss the implication of the ob- 
served faint end slope of the central galaxy luminos- 
ity and stellar mass functions for the relationship be- 
tween galaxies and dark matter halos. In Yang et al. 
(2005b) and Paper II, we have measured the CLFs for 
relatively massive halos with Mh > 10 12 ft~ 1 M Q directly 
from group samples. Unfortunately, galaxies in less mas- 
sive halos are not well studied, partly because the halo 
masses for small groups are difficult to estimate. Here, 
with the help of the L c — Mh (M* jC — Mh) relations that 
we obtained for relatively massive halos and using the lu- 
minosity functions and stellar mass functions for central 
galaxies, we can probe the halo properties of low-mass 
central galaxies in a statistical way. As discussed in Sec- 
tions 2] and [5l a small halo is usually dominated by a 
single central galaxy. The L c — Mh and M* jC — Mh re- 
lations for groups with masses down to ~ 10 118 ft _1 M Q 
can be obtained directly from the group catalogue, and 
are shown as open circles with error-bars in the upper 
panels of Fig. [5] For these relatively massive groups, we 
can also estimate cr(logL c ) and cr(logM», c ) as functions 
of halo mass. These are shown in the middle panels of 
Fig. [5]as the open circles with error-bars. Finally, the lu- 
minosity and stellar mass functions for central galaxies, 
which have been obtained in Section [3l are repeated in 
the lower panels of Fig. [H]as open circles with error-bars. 

Following paper II, we model the mean L c — Mh and 
M* iC — Mh relations using the following functional forms: 



and 



M* 



Mo 



(Mfe/MQ 
(1+Mfc/Mi)/ 3 ' 

(Mh/M^+P 



(19) 



(20) 



•{1 + Mh/MJP ' 
where Mi is a characteristic halo mass so that L c cx 
M^ +/3 (M*, c cx M£ +/3 ) for M h < Mi and L c cx M« 
(M* >c cx M%) for M h > Mi. Note that the parameters 
Mi, a and (3 may have different values in the two rela- 
tions. Clearly, for bright (massive) central galaxies, these 
relations are well constrained by the data points shown 
in the upper panels of Fig. [8l Unfortunately, no such 
direct constraints are available for the faint (low-mass) 
centrals. However, an indirect constraint comes from the 
observed luminosity and stellar mass functions of central 
galaxies. In general, we can write 



/■oo 

$(L C ) = / $ c en(L|Mfc 
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Fig. 8. — The properties of central galaxies: shown in the left and right panels are results for the luminosity L c and stellar mass M*, c , 
respectively. The open circles in each panel are the data we extracted from the SDSS DR4 group catalogue, while the solid, dotted and 
dashed lines are the best fitting results (see text for details). Upper- left panel: the L c — M h relation, where the data points are obtained 
from the left panel of Fig. 6 in Paper II. The long-dashed line in this panel illustrates the L c — Mf, relation predicted by Cacciato et al. 
(2008). Upper-right panel: the M», c — Mh relation, where the data points are obtained from the right panel of Fig. 6 in Paper II. Middle-left 
panel: the <r(logL c ) — Mh relation, where the data points are obtained from the lower-left panel of Fig. 4 in paper II. Middle-right panel: 
the <r(log Af*, c ) — M^ relation, where the data points are obtained from the lower-left panel of Fig. [5] Lower-left panel: the luminosity 
function of the central galaxies, $(L C ), where the data points are the same as those shown in the middle-left panel of Fig. [2] Lower-right 
panel: the stellar mass function of central galaxies, <J>(M» c ), where the data points are the same as those shown in the middle-right panel 
of Fig. [2] 



and 



*(M, iC ) 



<S> cen (M*\M h )n(M h )dM h , (22) 



where n(M) is the mass function of dark matter ha- 
los, and § ccn (L\M h ) and <& cen (M*\M h ) are the CLF and 
CSMF of central galaxies, respectively. If we model 
both Q ccn (L\Mh) and $> ccn (M*\Mh) with a lognormal 
form, they are completely determined by the mean re- 
lations, fH]) and (j2"0|) and the corresponding dispersions, 
er(logL c ) and cr(log Af». c ). As a simple model we first 
assume these dispersions to be constant: cr(logL c ) = (Tq 
and cr(logM» iC ) = cto (again Co may have different values 
in the two cases). Thus, for each case, we have five free 



parameters, Lq (or Mq), Mi, a, (3 and <jq. 

Using the halo mass function predicted by the ACDM 
'concordance' cosmology, we fit the models described 
above to the observational data shown in Fig. [U The 
fitting is performed with a standard least x 2 algorithm. 
For the luminosity of the central galaxies, we use y 2 — 
X 2 (L C ) + X 2 ( cr ) + X 2 (^(L C )), and for the stellar mass of 
the central galaxies, we use x 2 : 
2/r -\ Jit 



X 2 mM,, c )). Here, X 2 (L C ), X » 



X 2 (M» 



*,c) + X 2 {v) + 
and X 2 (HL C )) (or 

X 2 (M*, C ), X 2 ( a ) an d X 2 (^(M*,c))) are calculated accord- 
ing to the deviations of the model predictions from the 
observational data shown in the upper-, middle-, and 
lower-left (or right) panels of Fig. [51 respectively. Based 
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on all the data points shown in the left hand side panels 
of Fig. 11 we obtain logL = 9.9078, log Mi = 11.0096, 
a = 0.2566, /3 = 3.4037, cr = 0.1462. The resulting 
best fit is shown in each panel as the solid lines. The 
agreement between the best fit model and the data is 
remarkably good. The best fitting parameters indicate 
that L c oc M%- 7 for logM/j <C 11.0, suggesting that the 
star formation efficiency decreases dramatically in small 
halos. These results, especially the slopes of the L c - Mh 
relation (a ~ 0.25 at the massive end, and a + f3 » 1 
at the low-mass end), are in good agreement with pre- 
vious results (e.g., Vale & Ostriker 2004, 2006; Cooray 
2005; Yang et al. 2003; Yang et al. 2005c; van den Bosch 
et al. 2007; Hansen et al. 2007; Zheng et al. 2007; 
Popesso et al. 2007; Brough et al. 2008; Conroy & Wech- 
sler 2008). The physical reason for this change in slope 
is probably a combination of AGN feedback, and the 
change in the efficiency of radiative cooling, supernova 
feedback and dynamical friction (e.g. Lin et al. 2004; 
Cooray & Milosavljevic 2005). For comparison, we also 
show in the upper-left panel the mean L c — Mh relation 
obtained by Cacciato et al. (2008) based on the CLF 
models. The excellent agreement with our prediction in- 
dicates that although obtained via different methods, the 
mean L c — Mh relation are well constrained in both in- 
vestigations. 

For the stellar mass of the central galaxies, fitting the 
model to the data shown in the right panels of Fig. [5] 
gives log M = 10.3061, log Ml = 11.0404, a = 0.3146, 
(3 = 4.5427, cr = 0.1730. Note that M*, c oc M^ 9 for 
logM/j <C 11.0, which is even steeper than the L c - Mh 
relation. 

The above results are obtained under the assumption 
that the value of the dispersion, a, is independent of 
halo mass. Although consistent with the data for Mh > 
10 12 h~ 2 Mq, and supported by satellite kinematics and 
semi-analytical models (More et al. 2008), it is not well 
constrained in low mass halos. We therefore now test 
the impact of this assumption on our results. For this 
purpose, we consider two models where a is required to 
change from the observed value at log(M/j/ /i~ 1 M ) ~ 12 
either to or to 0.8 at log(M h / /i _1 M Q ) ~ 10 (see the 
dotted and dashed curves in the middle panels of Fig. 
[HI respectively). The best fits of these two test models 
are shown as the dotted and dashed lines in the other 
panels. Clearly, our results for the L c - Mh and M c - Mh 
relations are robust with respect to our assumption that 
a is independent of halo mass. 

7. SUMMARY 

In this paper, we have derived the luminosity and stel- 
lar mass functions for different populations of galaxies 
(central versus satellite; red versus blue; and galaxies in 
halos of different masses) , and for groups themselves, us- 
ing a large galaxy group catalogue constructed from the 
SDSS Data Release 4 (DR4). Our main results can be 
summarized as follows: 

1. For central galaxies, the conditional stellar mass 
function (CSMF), which describes the stellar mass 
distribution of galaxies in halos of a given mass 



can be well described by a log-normal distribution, 
with a width <7io g M, ~ 0.17, quite independent of 
the host halo mass. The median central stellar 
mass increases rapidly with halo mass, M» ex M^ 9 , 
for halos with masses Mh <C 10 11 H~ 1 Mq, but 
only slowly, M* oc Mfc 3 , for halos with M h > 
lO 13 ft _1 M . 

2. For satellite galaxies, the conditional stellar mass 
function in halos of different masses can be de- 
scribed reasonably well by a modified Schechter 
form than breaks away faster than the Schechter 
function at the massive end. The faint end slope 
appears to be steeper for more massive halos. On 
average, there are about 3 times as many central 
galaxies as satellites. 

3. When stellar mass functions are measured sepa- 
rately for galaxies of different colors, we find that 
the central population is dominated by red galaxies 
at the massive end, and by blue galaxies at the low- 
mass end. Among the satellite population, there 
are in general more red galaxies than blue ones. At 
the very low-mass end (M* < 10 9 /i~ 2 M Q ), there 
is a marked increase in the number of red centrals. 
We speculate that these galaxies are located close 
to large halos so that their star formation has been 
affected by their environments. 

4. The stellar-mass function of galaxy groups, which 
describes the number density of galaxy groups 
as a function of the total stellar mass of group 
member galaxies, is well described by a double 
power law, with a characteristic stellar mass at 
~ 4 x 10 10 /i~ 2 M©. This form is very different from 
that of the halo mass function, indicating that the 
efficiencies of star formation in halos of different 
masses are very different. 

5. The stellar mass function for the central galaxies 
can be used to provide stringent constraint on the 
mean M* jC - Mh relation for low-mass halos. 

We anticipate that a comparison of these results 
with predictions of numerical simulations and/or semi- 
analytical models will provide stringent constraints on 
how galaxies form and evolve in dark matter halos. 

ACKNOWLEDGMENTS 

We are grateful to the anonymous referee for useful and 
insightful comments that greatly helped to improve the 
presentation of this paper. XY acknowledges the Uni- 
versity of Massachusetts for hospitality where this work 
was finalized. This work is supported by the One Hun- 
dred Talents project, Shanghai Pujiang Program (No. 
07pjl4102), 973 Program (No. 2007CB815402), the CAS 
Knowledge Innovation Program (Grant No. KJCX2- 
YW-T05) and grants from NSFC (Nos. 10533030, 
10673023, 10821302). HJM would like to acknowledge 
the support of NSF AST-0607535, NASA AISR-126270 
and NSF IIS-0611948. 



REFERENCES 



Adelman-McCarthy J.K., et al., 2006, ApJS, 162, 38 



Baldry I.K., Glazebrook K., Brinkmann J., Ivezic Z., Lupton R.H., 
Nichol R.C., Szalay A.S., 2004, ApJ, 600, 681 



18 



Yang et al. 



Baldry I.K., Glazebrook K., Driver S.P., 2008, MNRAS, 388, 945 
Bell E.F., Mcintosh D.H., Katz N., Weinberg M.D., 2003, ApJS, 
149, 289 

Berlind A.A., Weinberg D.H., 2002, ApJ, 575, 587 
Berlind A.A., et al. 2006, ApJS, 167, 1 

Berlind A.A., Kazin E., Blanton M.R.. Pueblas S., Scoccimarro R., 

Hogg D.H., 2007, preprint ( astro-ph/0610524 ) 
Blanton M. R. et al. , 2003a, AJ, 125, 2348 
Blanton M. R. et al. , 2003b, ApJ, 592, 819 

Blanton M.R., Eisenstein D.J., Hogg D.W., Schlegel D.J., 

Brinkmann J., 2005a, ApJ, 629, 143 
Blanton M.R. et al. , 2005b, AJ, 129, 2562 
Blanton M.R., Roweis S., 2007, AJ, 133, 734 
Borch A., et al., 2006, A&A, 453, 869 

Brough S., Couch W.J., Collins C.A., Jarrett T., Burke D.J., Mann 

R.G., 2008, MNRAS, 385L, 103 
Brown M.J. I., et al., 2008, ApJ, 682, 937 

Bullock J.S., Wechsler, R.H., Somerville R.S., 2002, MNRAS, 329, 
246 

Cacciato M., van den Bosch F.C ., More S., Li R., Mo H.J., Yang 
X., 2008, MNRAS, submitted (arXiv:0807.4932 ) 

Coil A.L., et al., 2006, ApJ, 638, 668 

Colless M., et al., 2001, MNRAS, 328, 1039 

Collister A.A., Lahav O., 2005, MNRAS, 361, 415 

Conroy C, Wechsler R.H., 2008, preprint, arXiv:0805.3346 

Cooray A., 2005, MNRAS, 364, 303 

Cooray A., Milosavljevic M., 2005, ApJ, 627, L89 

Cooray A., 2006, MNRAS, 365, 842 

Fontana A., et al., 2006, AfeA 459, 745 

Guo Y., et al., 2009, preprint, arXiv:0901.1150 

Hansen S.M., Sheldon E.S., Wechsler R.H., Koester B.P., 2007, 
preprint ( arXiv:0710.3780 ) 

Jing Y.P., Mo H.J., Borner C, 1998, ApJ, 494, 1 

Jing Y.P., Borner C, Suto Y., 2002, ApJ, 564, 15 

Kassin S.A., et al., 2007, ApJ, 660L, 35 

Kroupa P., 2001, MNRAS, 322, 231 

Li C, Kauffmann C, Jing Y.P., White S.D.M., Borner C, Cheng 

F.Z., 2006, MNRAS, 368, 21 
Lin H., Kirshner R.P., Shectman S.A., Landy S.D., Oemler A., 

Tucker D.L., Schechter P.L., 1996, ApJ, 464, 60 
Lin Y.T., Mohr J. J., 2004, ApJ, 617, 879 

Ludlow A.D., Navarro J.F., Springel V., Jenkins A., Frenk C.S., 

Helmi A., 2008, preprint (arXiv:0801.1127) 
Madgwick D.S., et al., 2002, MNRAS, 333, 133 
McGaugh S.S., 2005, ApJ, 632, 859 

Mcintosh D.H., Guo Y., Hertzberg J., Kat z N.. Mo H.J van den 
Bosch F.C., Yang X., 2007, preprint (arXiv:0710.2157) 

Mo H.J., Yang X.H., van den Bosch F.C., Jing Y.P., 2004, MNRAS, 
349, 205 

More S. van den Bosch F.C., Cacciato M., Mo H.J., Yang X., Li 

R.. 2008, MNRAS, submitted (arXiv:0807:4532) 
Norberg P., et al., 2002, MNRAS, 336, 907 

Panter B., Jimenez R., Heavens A. F., Chariot S., 2007, MNRAS, 
378, 1550 

Peacock J. A., Smith R.E., 2000, MNRAS, 318, 1144 
Petrosian V., 1976, ApJ, 209, LI 
Pizagno J. et al., 2007, AJ, 134, 945 

Popesso P., Biviano A., Bohringer H., Romaniello M., 2006, A&A, 
445, 29 

Robotham A., Wallace C, Phillipps S., De Propris R., 2006, ApJ, 
652, 1077 

Schlegel D.J., Finkbeiner D.P., Davis M., 1998, ApJ, 500, 525 

Scoccimarro R., Sheth R.K., Hui L., Jain B., 2001, ApJ, 546, 20 

Scranton R., 2002a, MNRAS, 332, 697 

Seljak U.,2000, MNRAS, 318, 203 

Sheldon E.S., et al., 2007, preprint, arXiv:0709.1162 

Skibba R.A., Sheth R.K., Martino M.C., 2007, MNRAS, 382, 1940 

Skibba R.A., Sheth R.K., 2008, preprint, arXiv:0805.0310 

Skibba R.A., 2008, preprint, arXiv:0805.1233 



Spergel, D. N., et al. 2007, ApJS, 170, 377 

Stoughton C, et al., 2002, AJ, 123, 485 

Strauss M.A., et al., 2002, AJ, 124, 1810 

Tago E., et al., 2006, Astronomische Nachrichten, 327, 365 

Tempel E., Einasto J., Einasto M., Saar E., 2008, preprint, 

[arXiv:0805.4264] 
Tinker J.L., Weinberg D.H., Zheng Z., Zehavi I., 2005, ApJ, 631, 

41 

Tully R. B., Pierce M.J., 2000, ApJ, 533, 744 

Vale A., Ostriker J. P., 2004, MNRAS, 353, 189 

Vale A., Ostriker J. P., 2006, MNRAS, 371. 1173 

Vale A., Ostriker J. P., 2007, preprint (astro-ph/0701096 ) 

van den Bosch F.C., Yang X., Mo H.J., 2003, MNitAb, 340, 771 

van den Bosch F.C., Mo H.J., Yang X., 2003, MNRAS, 345, 923 

van den Bosch F.C., Weinmann S.M., Yang X., Mo H.J., Li O, 

Jing Y.P., 2005, MNRAS, 361, 1203 
von der Linden A., Best P.N., Kauffmann C, White S.D.M., 2007, 

MNRAS, 379, 867 

van den Bosch F.C., Yang X., Mo H.J., Weinmann S.M., Maccio 
A., More S., Cacciato M., Skibba R., Kang X., 2007, MNRAS, 
376, 841 

van den Bosch F.C., Aquino D., Yang X., Mo H.J., Pasquali A., 
Mcintosh D.H., Weinmann S.M., Kang X., 2008, MNRAS, 387, 
79 

Wang Y., Yang X., Mo H.J., van den Bosch F.C., Chu Y., 2004, 

MNRAS, 353, 287 
Wang Y., Yang X., Mo H.J., van den Bosch F.C., Weinmann S.M., 

Chu Y., 2008a, ApJ, 687, 919 
Wang Y., Yang X., Mo H.J., van den Bosch F.C., Katz N, 

Pasq uali A.. Mc intosh D.H., Weinmann S.M., 2008b, preprint, 

[aFXiv:0812.3723] 
Warren M.S., Abazajian K., Holz D.E., Teodoro L., 2006, ApJ, 

646, 881 

Weinmann S.M., van den Bosch F.C., Yang X., Mo H.J., 2006a. 
MNRAS, 366, 2 

Weinmann S.M., van den Bosch F.C., Yang X., Mo H.J., Croton 

D.J., Moore, B., 2006b, MNRAS, 372, 1161 
White M., Zheng Z., Brown M.J.I., Dey A., Jannuzi B.T.. 2007. 

ApJ, 655, 69 

Yan R., Madgwick D.S., White M., 2003, ApJ, 598, 848 

Yan R., White M., Coil A.L., 2004, ApJ, 607, 739 

Yang X., Mo H.J., van den Bosch F.C., 2003, MNRAS, 339, 1057 

Yang X., Mo H.J., Jing Y.P., van den Bosch F.C., Chu Y.Q., 2004, 

MNRAS, 350, 1153 
Yang X., Mo H.J., van den Bosch F.C., Jing Y.P., 2005a, MNRAS, 

356, 1293 

Yang X., Mo H.J., van den Bosch F.C., Jing Y.P., 2005b, MNRAS, 

357, 608 

Yang X., Mo H.J., Jing Y.P., van den Bosch F.C. 2005c, MNRAS, 

358, 217 (Y05c) 

Yang X., Mo H.J., van den Bosch F.C., Weinmann S.M., Li O, 

Jing Y.P., 2005d, MNRAS, 362, 711 
Yang X., Mo H.J., van den Bosch F.C., 2006, ApJ, 638L, 55 
Yang X., Mo H.J., van den Bosch F.C., Pasquali A., Li O, Barden 

M., 2007, ApJ, 671, 153 (Paper I) 
Yang X., Mo H.J., van den Bosch F.C., 2008a, ApJ, 676, 248 (Paper 

II) 

Y ang X., Mo H .J., van den Bosch F.C., 2008b, preprint, 
arXiv:0808.2526 

Zandivarez A., Martnez H.J., Merchan M.E., 2006, ApJ, 650, 137 

Zehavi I., et al., 2002, ApJ, 571, 172 

Zehavi I., et al., 2004 ApJ, 608, 16 

Zehavi I., et al., 2005 ApJ, 630, 1 

Zheng Z., et al., 2005, ApJ, 633, 791 

Zheng Z., Coil A.L., Zehavi I., 2007, ApJ, 667, 760 

Zheng Z., Weinberg D.H., 2007, ApJ, 659, 1 



